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DYNAMIC  ICE-STRUCTURE  INTERACTION 
DURING  CONTINUOUS  CRUSHING 


Mauri  Maattanen  '■ 
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INTRODUCTION 

The  interactions  between  the  ice  and  the  struc¬ 
ture  when  a  constant-thickness  ice  sheet  crushes 
against  a  vertical  pile  have  long  been  studied,  both 
in  the  laboratory  and  in  the  field.  However,  due 
to  the  limits  of  the  test  arrangements,  the  great 
number  of  parameters  to  be  considered,  and  the 
complexity  of  the  problem,  dynamic  ice-structure 
interactions  still  are  not  fully  understood. 

Most  of  the  research  in  this  field f  - 1 0  has  focused 
on  the  maximum  ice  force;  dynamic  variations  have 
received  little  or  no  attention.  Most  researchers  try 
to  avoid  dynamic  ice-structure  interactions  by  mak¬ 
ing  the  loading  apparatus  as  stiff  as  possible.  In 
addition  the  stroke  of  the  load  actuator  has  often 
been  so  short  that  continuous  crushing  interaction 
is  impossible.  Only  one  experiment1 1  has  tested 
dynamic  ice-structure  interactions. 

In  laboratory  tests  it  is  easy  to  control  and  vary 
parameters  and  to  repeat  the  test  if  necessary.  The 
difficulty  in  laboratory  tests  is  scaling  the  properties 
of  both  the  ice  and  the  structure  correctly.  It  is  not 
possible  to  scale  down  all  ice  properties  simultan¬ 
eously.1 1  The  effects  of  scaling  on  the  relationship 
between  ice  strength  and  loading  rate  are  unknown. 

Field  tests  give  valuable  data  on  real  ice-structure 
interaction  without  any  scale  effects.  The  interpre¬ 
tation  of  the  results,  however,  is  more  difficult. 
Usually  it  is  not  possible  to  measure  the  basic  param¬ 
eters  of  ice  (velocity,  thickness,  temperature,  and 
grain  size  and  orientation)  while  the  ice  is  moving. 
Published  reports5'8'  12'14only  give  estimates  of 
these  parameters.  The  chances  of  changing  one 
parameter  in  the  field  are  minimal. 


There  are  different  opinions  about  the  participa¬ 
tion  of  the  structure  itself  in  the  interaction:  most 
of  the  authors7-18'14'16  explained  that  the  dynamic 
ice  force  results  from  the  properties  of  ice  alone, 
while  Maattanen15  proposed  that  the  dynamic  re¬ 
sponse  of  the  structure  affects  the  resulting  ice  force 
and  ice  failure  process,  so  that  the  properties  of  both 
the  ice  and  the  structures  are  significant.  With  very 
stiff  structures  the  effect  of  structural  response  be¬ 
comes  insignificant,  but  with  slender  pile  and  pier 
structures  the  dynamic  stiffness  plays  an  important 
role  in  the  resulting  interacting  ice  force.  Earlier 
tests  lacked  the  power  of  modern  signal-analyzing 
techniques  and  the  versatility  of  computer  modeling. 
Now  the  dynamic  properties  of  the  structures  related 
to  primary  structural  parameters  (the  width  and 
shape  of  the  area  exposed  to  ice  actions)  can  also  be 
assessed. 

There  have  been  only  two  theoretical  approaches 
to  the  ice-structure  interaction  problem.  Reddy  ct 
al.17  used  a  statistical  approach  in  which  the  ice 
force  excitation  (assumed  to  be  random)  was  mod¬ 
ified  by  the  properties  of  the  structure  to  give  the 
dynamic  response  spectra.  In  this  forced  vibration 
approach  no  real  ice-structure  interaction  occurs. 
Maattanen15  related  the  ice  strength  dependence  on 
the  loading  rate  to  the  deformation  rate  of  the  struc¬ 
ture  at  the  ice  action  point.  Hence,  the  ice  and  the 
structure  form  a  closed  loop  system:  each  affects 
the  other.  With  this  interaction  model  it  is  possible 
to  solve  the  resulting  ice  force  vs  time  function  with¬ 
out  any  arbitrary  assumptions,  because  only  known 
physical  properties  of  the  ice  and  the  structure  are 
used.  The  least  known  parameter  is  the  dependence 
of  the  ice  strength  on  the  loading  rate.  There  exist 


TEST  ARRANGEMENTS 


only  a  few  measurements  of  this  dependence.7*10*14 
Ice  strength  is  a  statistical  quantity,  and  hence  the 
published  curves  of  crushing  strength  vs  loading  rate 
are  averages. 

This  research  is  a  continuation  of  the  dynamic  ice- 
structure  interaction  tests  that  started  in  1978  at  the 
CRREL  Ice  Engineering  Facility.  In  its  present  form 
the  laboratory  arrangements  allow  a  number  of 
parameters  to  be  measured:  ice  and  structural  prop¬ 
erties  in  a  scale  from  1 :4  to  1 :40,  temperature  ranges 
from  0°  to  -20°C,  and  ice  velocities  that  yield  load¬ 
ing  rates  for  ice  failure  from  ductile  to  brittle.  The 
first  report11  describes  dynamic  ice-structure  inter¬ 
action  test  arrangements  and  procedures,  instrumen¬ 
tation,  modeling  laws  and  analysis  methods.  Test 
results  at  high  velocity  range  were  given.  At  the 
most  interesting  low-velocity  range  (the  transition 
from  ductile  to  brittle  ice  failure  mode),  results  were 
unreliable  due  to  jerking  ice  movement,  which  was 
caused  by  an  ice-pulling  mechanism  that  was  too 
flexible. 

The  main  objective  of  this  research  was  to  meas¬ 
ure  the  dependence  of  ice  forces  and  crushing  fre¬ 
quency  on  the  properties  of  both  the  ice  and  the 
structure  during  continuous  crushing  at  the  low  ve¬ 
locity  range.  The  varied  parameters  include:  ice 
thickness,  pile  diameter  at  the  waterline,  pile  natural 
frequencies  and  natural  modes.  The  ice-induced 
positive  and  negative  damping  effects  and  the  stabil¬ 
ity  and  limit  cycles  of  self-excited  vibrations  were 
also  emphasized  in  the  data  analysis. 


The  simulation  of  moving  ice  crushing  against  a 
bottom-founded  pile  was  arranged  in  the  CRREL  Ice 
Engineering  Facility  so  that  the  scale  model  of  the 
structure  was  fixed  to  the  bottom  of  the  test  basin 
and  the  ice  was  pushed  by  a  carriage  bridge.  The 
velocity  of  the  ice  was  controlled  by  a  closed  loop 
control  system.  The  normal  procedure  was  to  main¬ 
tain  a  constant  acceleration  in  order  to  sweep  the 
whole  velocity  range  in  one  run. 

The  actual  structures  could  be  scaled  with  an  ex¬ 
act  geometrical  similitude.  Then,  however,  the 
structural  properties  would  be  difficult  to  change  in 
subsequent  tests.  For  studying  ice-structure  inter¬ 
actions  it  is  enough  to  scale  down  correctly  only 
those  natural  modes  that  contribute  significantly  to 
deformations  at  the  ice  action  point.  Therefore,  a 
test  pile  and  supporting  structure  (Figs.  1  and  2)  was 
designed  so  that  the  stiffness,  natural  frequencies 
and  modes  could  be  easily  varied.  The  width  at  the 
waterline  and  the  shape  of  the  area  affected  by  the 
ice  could  be  varied  by  changing  the  collar  at  the 
waterline.  By  this  method  most  of  the  single-pile, 
bottom-founded  structures  could  be  simulated. 

The  different  configurations  of  the  test  pile  were 
calibrated  before  testing.  The  transfer  function 
measurement  gives  all  the  information  about  the 
dynamic  properties  of  the  structure.  This  calibration 
was  carried  out  by  impacting  the  pile  at  the  water¬ 
line  with  a  load  cell.  From  the  ratio  between  the 
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Figure  2.  Test  pile  blueprint. 


output  signal  at  the  strain  gauges  in  the  test  pile 
and  the  input  impact  load  signal,  the  transfer  func¬ 
tion  was  calculated.  Natural  frequencies  and  modal 
damping  are  then  readily  available,  and  the  static 
stiffness  against  the  ice  force  can  be  calculated. 
Figure  3  shows  the  transfer  and  coherence  functions 
for  Tests  34-51 .  In  addition,  damping  was  measured 
from  the  decay  rates  of  the  acceleration  signals  after 
impacts.  The  ranges  of  the  available  first  and  second 
natural  frequencies  without  and  with  a  superstruc¬ 
ture  tuning  weight  are  given  in  Figure  4.  Also,  the 
curve  for  static  stiffness  is  presented  as  a  function 
of  test  pile  support  location. 

The  instrumentation  in  the  test  pile  included 
strain  gauges  used  to  measure  the  reaction  of  the 
ice  force  at  the  upper  support  clamp.  In  the  super¬ 
structure  just  above  the  waterline  were  two  sets  of 
strain  gauges  used  to  measure  the  bending  moment 
and  shear  force,  both  of  which  were  used  to  elimin¬ 
ate  the  inertial  effects  of  the  superstructure  on  the 
ice  force.  Accelerometers  on  the  superstructures 
could  be  attached  at  different  levels  and  in  different 


directions;  they  were  used  to  measure  the  dynamic 
response.  In  some  tests  one  accelerometer  was  in¬ 
stalled  on  the  ice  or  on  the  carriage  for  measuring 
ice  jerking.  The  instruments  for  recording  and  an¬ 
alyzing  data  and  the  analyzing  procedures  are  the 
same  as  explained  by  Maattanen. 11  Figure  5  shows 
the  signal  recording  and  processing  block  diagram. 

As  an  improvement  to  the  first  test  setup,11  the 
test  carriage  bridge  was  used  to  pull  the  ice  in  the 
test  basin.  However,  neither  of  its  test  setups  was 
capable  of  providing  velocities  low  enough  for  ice- 
structure  interaction  tests.  Therefore,  the  main 
drive  of  the  test  carriage  was  disconnected  and  its 
reduction  gear  removed;  they  were  replaced  by  the 
DC-drive  and  reduction  gear  from  the  previous  test 
setup.11  Then  the  reduction  gear  was  driving  the 
drive  shafts  of  the  carriage  instead  of  winches.  This 
gave  a  controllable  speed  range  from  4  to  104  mm/s. 
With  the  scaling  factors  the  loading  rate  correspond¬ 
ing  to  the  transition  in  ice  strength  from  ductile  to 
brittle  was  about  10  mm/s.  The  average  stress  rate 
that  could  be  achieved  ranged  from  80  to  6000  kPa/s. 
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Clomp  Hole 

Figure  4.  Test  pile  natural  frequencies  and  stiffness 
vs  support  damp  attachment  position. 

The  only  shortcoming  of  the  test  arrangements 
was  the  stiffness  of  the  ice-pulling  system,  despite  the 
improvement  over  the  previous  system.  It  was  pos¬ 
sible  to  carry  out  tests  throughout  the  velocity  range. 
Even  though  the  stiffness  of  the  bridge  itself  was 
adequate,  the  torsional  rigidity  of  the  drive  shafts 
reduced  the  theoret:cal  values  of  the  lowest  natural 
frequency  of  the  bridge  to  below  5.6  Hz  and  reduced 
the  measured  value  to  2.5  Hz.  Unfortunately  this 
frequency  is  in  the  range  of  interest  for  the  test  pile 
and  interferes  with  the  resulting  ice  force.  With  the 
low  velocities  in  some  tests  the  ice  velocity  was  not 
constant  but  included  a  carriage  vibration  compo¬ 
nent.  For  future  tests  it  may  be  better  to  make  the 
test  pile  more  flexible  and  the  ice  weaker  to  minimize 
the  effect  of  carriage  flexibility  on  the  ice  force. 


ICE  PROPERTIES 

Scaling  down  the  ice  properties  required  that  all 
the  parameters  affecting  the  studied  phenomenon 
be  correctly  taken  into  account.  During  studies  of 
ice  crushing  against  a  pile,  geometric  similitude  of 
the  pile  diameter,  the  thickness,  strength,  stiffness 
and  failure  patterns  of  the  ice,  and  the  grain  size  and 
structure  are  necessary. 

Six  ice  fields  were  frozen  and  tested  in  the  test 
basin  during  December  1979  and  January  1980.  The 
first  was  just  for  cleaning  the  dust  and  dirt  from  the 
water  surface  of  the  tank.  The  instrumentation  was 
tested  as  this  ice  sheet  was  pulled  across  the  basin. 
The  second  ice  field  was  made  of  freshwater  ice. 

The  seeding  was  done  with  a  water  spray;  the  results 
were  poor,  with  several  patches  of  ice  having  uncon¬ 
trolled  large  crystals.  The  freshwater  results  arc 
therefore  not  reliable. 

The  use  of  urea  instead  of  natrium  chloride  to 
scale  down  ice  properties  has  been  shown  to  be 
promising  in  laboratory  tests.18  To  avoid  the  severe 
corrosion  problems  with  salt,  urea  was  used.  For 
about  a  1:10  reduction  in  ice  strength  a  0.6%  urea 
concentration  was  chosen  for  use  in  the  next  four 
ice  fields.  This  appeared  to  be  a  good  choice. 

To  produce  the  urea-doped  ice  the  seeding  was 
done  using  steam  at  an  air  temperature  of  -1  5°C, 
yielding  fine  initial  grains.  About  15  mm  from  the 
top  a  columnar-grained  ice  with  a  horizontal  c-axis 
started  to  grow.  As  crystallographic  studies  show 
(Fig.  6),  the  urea  ice  is  similar  to  ser  ice,  fulfilling 
the  geometric  similitude  requirement.  The  cracking 
patterns  (radial  cracks  both  longitudinal  and  20-40° 
off  the  centerline)  were  typically  similar  to  those  of 


Figure  5.  Signal  recording  and  processing  equipment. 
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Figure  6.  Vertical  cross  section  (rectangle)  and  horizontal  cross  sections  (circles)  of  0.6 % 
urea  ice  ( tests  3  7-40).  The  grid  size  is  5x5  mm. 


Figure  7.  Crushing  truth  trass  sections. 
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saline  ice.  Also,  the  cross  sections  of  the  crushing 
track  (Fig.  7)  were  similar  to  those  of  saline  ice. 
However,  it  appears  that  ice  thicknesses  well  above 
20  mm  shouid  be  used  so  that  the  proportion  of 
columnar  crystals  will  be  sufficient.  The  strength 
of  urea  ice  reduces  during  the  warm-up  period, 
which  has  to  be  longer  than  with  saline  ice. 

For  crushing,  the  best  measure  of  strength  would 
be  plane  compressive  or  shear  strength.  Due  to  the 
lack  of  procedures  for  testing  these  properties, 
standard  cantilever  beam  flexural  strength  and  plate 
bending  tests  were  used  as  index  tests.21  The  ice 
was  warmed  by  keeping  the  air  temperature  at  -1°C 
until  the  required  flexural  strength  was  achieved. 
Then  the  crushing  test  with  a  model  structure  was 
started;  the  crushing  strength  was  calculated  from 
the  measurement  data.  The  warm-up  period  from 
an  ambient  temperature  of  -15°C  usually  took  about 
six  hours.  The  C/a.fmodulus  of  elasticity/flexural 
strength)  ratio  remained  well  about  2000.  The  in¬ 
finite  plate  bending  modulus  of  elasticity  tests  were 
in  the  process  of  being  developed;  the  displacement 
transducers  were  located  too  far  from  the  load 
actuator.  Therefore,  the  modulus  of  elasticity  val¬ 
ues  are  not  very  accurate,  and  there  were  no  meas¬ 
urements  at  all  for  the  last  two  ice  fields. 


Table  1  gives  the  main  properties  of  the  ice  fields 
in  these  tests  Mote  that  the  thinner  ice  was  stronger 
than  the  thicker  ice.  Ice  sheet  1 5  was  the  weakest 
and  ice  sheet  16  was  the  strongest  in  crushing,  even 
though  the  differences  in  bending  strength  were 
small.  This  occurs  with  urea  ice  because  the  top 
layer  of  fine,  randomly  oriented  ice  crystals  domi¬ 
nates  flexural  strength,  while  the  whole  thickness  is 
as  important  in  crushing.  During  warm-up  the 
strength  in  the  columnar  section  deteriorates  faster 
than  in  the  randomly  oriented  section.  The  columns 
start  to  get  loose  from  each  other,  and  the  ice  crushes 
easier.  This  phenomenon  is  also  observed  in  the 
laboratory  with  sea  ice. 

CRUSHING  PATTERNS 

Dynamic  ice-structure  interactions  appeared  in 
each  of  the  tests.  Especially  good  repetitive  ice 
force  fluctuations  were  achieved  with  urea  ice.  The 
freshwater  ice  exhibited  more  randomness,  evidently 
Uue  to  its  more  brittle  nature.  As  the  velocity  of  ice 
increased,  the  force  fluctuation  pattern  changed  sig¬ 
nificantly. 
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Table  1 .  Ice  properties  in  tests  32-51 . 

Pile  Ice  Ice  Flexural  Modulus  Crushing 

Ice  no.  Test  diameter  thickness  temperature  strength  of  elasticity  strength 
(%  urea)  no.  (mm) _ (mm)  (°C) _ (kPg) _ (GPa) _ (MPa) _ Notes 


11 

32 

102 

(0.0) 

33 

102 

12 

34 

102 

(0.0) 

35 

178 

36 

64 

37 

64 

13 

38 

102 

(0.60) 

39 

178 

40 

178 

41 

178 

14 

42 

102 

(0.60) 

43 

102 

44 

178 

45 

102 

15 

46 

178 

(0.60) 

47 

178 

48 

64 

16 

49 

64 

(0.60) 

50 

102 

51 

178 

30 

- 

30 

- 

35 

-0.2 

38 

-0.1 

830 

38 

-0.1 

50 

-1.0 

51 

-0.8 

150 

52 

-0.6 

53 

-0.4 

53 

-0.9 

53 

-0.8 

120 

54 

-0.8 

55 

-0.7 

79 

79 

79 

82 

160 

Layer  of  dust  above  water 


3.26  Partly  unsuccessful  seeding 

1.69  by  water  spray 

3.30 


0.94  Some  fluctuations  of  air 

3.90  0.80  temperature  during  freecing 

0.92  -20°C  ...  -10°C 

0.86 


0.97  Severe  fluctuations  of  air 

3.20  1.07  temperature  during  freezing 

1.36  -20°C  ...  -4°C 

1.02 


0.71 

0.71 

0.64 

0.69 


22 

21 

20 


150 


1.4  3 
1.48 

1.5  3 


At  low  speeds  a  sawtooth  pattern  formed.  The 
ice  edge  stuck  to  the  pile  and  deflected  it  with  al¬ 
most  the  velocity  of  the  ice  movement.  The  low 
loading  rate  initially  caused  ductile  failure  of  the 
ice,  but  the  spring-back  rate  of  the  structure  as  the 
ice  failed  by  crushing  made  the  failure  brittle.  Due 
to  the  inertia  of  the  structure  the  spring-back  ex¬ 
ceeded  the  equilibrium  position.  After  the  maximum 
spring-back  the  pile  started  to  spring  forward  in  the 
direction  of  ice  movement.  At  low  ice  velocities 
the  spring-forward  rate  was  higher  than  the  ice  ve¬ 
locity,  causing  the  structure  to  pull  away  from  the 
ice.  Shortly,  however,  the  edge  of  the  ice  hit  again, 
causing  another  cycle  of  crushing  and  separating. 
Because  the  ice  properties  were  random,  these  load¬ 
ing  cycles  varied;  however,  the  overall  patterns  re¬ 
main  the  same.  Good  examples  are  the  constant  ice 
velocity  test  40  (Fig.  8a)  and  a  sample  of  fest  44 
(Fig.  8b). 

With  the  increasing  ice  velocity  the  sawtooth 
pattern  persisted  but  included  more  randomness. 
After  the  spring-back  the  following  spring-forward 


rate  no  longer  exceeded  the  ice  velocity;  thus,  the 
contact  between  ice  and  structure  was  continuous, 
and  the  ice  force  no  longer  dropped  to  zero  after 
crushing.  During  the  ice  force  build-up  the  loading 
rate  can  still  cause  ice  failure  in  the  ductile  range, 
even  though  the  ice  velocity  alone  would  make  the 
failure  brittle.  The  deflection  rate  of  the  pile  affects 
the  relative  velocity  between  the  ice  and  pile,  and  it 
is  this  relative  velocity  that  determines  the  momen¬ 
tary  loading  rate. 

The  average  ice  velocity  has  its  effects,  of  course. 
The  maximum  ice  forces  usually  occurred  at  the  low 
ice  velocities  corresponding  to  the  ductile  loading 
rate.  In  some  instances  the  maximum  ice  forces 
appeared  at  the  high  velocities  corresponding  to  the 
brittle  loading  rate.  Whether  this  was  a  result  of 
randomness  in  ice  strength  and  contact  or  still  a  re¬ 
sult  of  low  relative  velocity  is  difficult  to  judge.  In 
most  cases  the  high  ice  force  at  the  beginning  of  the 
test  was  due  to  better  contact  between  the  ice  and 
the  pile  than  after  crushing  caused  random  contact 


areas. 


100  -  150mm 


Figure  9.  Crushing  tracks  and  cracking  patterns  from  above. 


The  failure  of  the  ice  before  the  circular  pile 
caused  radial  cracks  to  propagate.  Only  in  test  41 
were  circumferential  cracks  observed.  Longitudinal 
cracks  were  the  most  common,  but  cracks  aligned 
20-40°  aside  were  also  frequent.  The  lengths  of 
the  cracks  were  from  0.3  to  3  meters.  The  longest 
one  split  the  ice  field  into  two  parts.  In  a  few  cases 
transverse  cracks  formed  on  both  sides  of  the  test 
pile,  cutting  the  moving  ice  field.  The  cracks  appar¬ 
ently  had  little  cfTect  on  the  ice  crushing  itself.  The 
measurement  data  also  give  no  clear  indication  of 
the  reduction  in  ice  forces  due  to  cracks.  The  ob¬ 
vious  explanation  for  this  reduction  is  that  the  edges 
of  the  test  basin  conf  ined  further  lateral  movement 
of  the  ice,  and  thus  the  only  way  for  ice  to  pass  the 
pile  was  to  be  crushed. 

f  urther  details  of  the  ice  failure  process  in  crush¬ 
ing  were  studied  by  taking  a  16  mm  film.  The  re¬ 
sults  were  not  satisfactory,  and  more  information 
was  gained  by  observing  the  cross  sections  of  crush¬ 
ing  tracks  and  broken  pieces  of  ice  directly.  Figure 
7  shows  the  cross  sections  of  the  c  ushing  tracks; 
f  iguie  9  shows  some  top  views. 

The  most  common  crushing  track  for  freshwater 
ice  was  tapered.  The  taper  angle  was  not  constant 
but  varied  randomly  from  vertical  to  30°  and  on  a 
few  occasions  to  45°.  Wedged  cross  sections  also 
appeared  frequently.  The  width  of  the  track  was 
almost  exactly  the  same  as  the  width  of  the  pile. 

The  edges  at  the  narrowest  point  of  the  cross  sections 
were  smoothed  because  the  ice  yielded  under  pres¬ 
sure  aeainst  the  sides  ol  the  pile. 

I  he  crushing  tracks  of  urea  ice  were  mostly  sim¬ 
ilar  to  those  of  freshwater  ice.  A  different  feature 
was  a  lipped  track,  which  also  occurs  in  saline  ice. 

This  pattern  resulted  from  the  deterioration  of  bonds 
between  columnar  grains  during  the  waim-up  period 
before  the  tests.  During  crushing,  as  the  crushed  ice 
rubble  was  pushed  aside,  it  loosened  entire  columnar 
grains  instead  of  splitting  them  as  in  the  more  com¬ 
mon  tapered  pattern.  In  rare  cases  the  ice  split  at 
the  center.  If  the  horizontal  center  crack  developed 
in  the  zone  of  transition  from  randomly  oiiented 
icecivstals  tocoF'mnar  grains,  it  would  cieate  the 
lipped  ciushing  liack.  In  the  wedged  track,  wedges 


extending  up  to  50  mm  were  split  randomly  from 
both  upper  and  lower  surfaces. 

Top  views  of  the  crushing  patterns  are  shown  in 
Figure  9.  The  most  common  type  had  radial  cracks; 
however,  sometimes  no  radial  cracks  were  present 
and  in  some  cases  a  wedge  of  ice  in  front  of  the  pile 
crushed  into  small  pieces  at  once.  There  was  no 
correlation  between  the  cracking  patterns  and  ice 
velocity.  With  thin  ice,  wide  diameter  pile  and  low 
velocity,  ice  failed  by  buckling  (tests  50  and  51 ). 

The  pieces  of  crushed  ice  were  usually  very  small, 
typically  below  5  mm;  only  the  wedges  of  the  type 
shown  in  Figure  7  were  larger. 

All  of  these  crushing  and  cracking  patterns  except 
the  lipped  crushing  track  have  been  observed  in  the 
field  by  the  author.  This  confirms  the  usefulness 
of  laboratory  tests  in  simulating  ice  crushing  against 
piles. 


MAXIMUM  ICE  FORCE  VS  VELOCITY 

Table  2  is  a  compilation  of  the  maximum  iccordcd 
ice  forces  and  their  corresponding  velocities.  The 
ice  lorccs  are  refined  from  the  dynamic  inertial  ef¬ 
fects  of  the  superstructure  by  using  the  bending 
moment  and  shear  force  measurements  above  the 
waterline  in  the  test  pile.  An  error,  smaller  than  that 
ol  the  superstructure  elfects,  remains  and  is  attrib¬ 
uted  to  the  inertial  elfects  ol  the  structures  below 
the  waterline.  In  tests  34,  35  and  36,  a  floppy-disk 
memory  error  makes  the  data  unreliable  for  velocities 
over  20mm/s. 

The  data  show  that  hall  ol  the  maximum  occurred 
in  tile  low  velocity  range  and  half  in  the  high  velocity 
range.  However,  the  overall  number  ol  occurrences 
in  each  test  in  the  low  velocity  lange  is  more  fiequent. 
At  high  velocities  the  ice  lorce  is  more  landom,  and 
only  when  a  combination  of  coned  dynamic  ampli¬ 
fication  rate,  good  ice  contact  and  high  icestienglh 
occur  will  the  resulting  ice  force  exceed  the  ice  loices 
rn  the  low  velocity  range.  Considering  the  average 
of  averaged  ice  forces,  then'  is  a  reduction  of  38—15'.'. 
fiom  low  to  high  veloc  its  and  a  reduction  of  I  I  - 1 4% 
loi  average  ol  maximums  m  tests  37-15  I  m  tests 
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Table  2.  Maximum  ice  force  vs  velocity. 


Test 

number 

F  * 

max 

(kN) 

vot  ^  max 
(mm/s) 

Fav  max 
(kN) 

vat  Fov  max 
(mm/s) 

Fsrtf 

(kN) 

34 

19.3 

18 

8.2 

17 

5.7 

35 

16.2 

108 

2.6 

14 

1.9 

36 

15.7 

102 

2.0 

18 

0.9 

37 

6.6 

11 

5.3 

9 

0.9 

38 

9.1 

30 

4.9 

14 

1.3 

39 

14.5 

11 

7.1 

9 

4.3 

40 

14.4 

12 

9.6 

6 

2.4 

41 

14.5 

1 6 

9.1 

18 

3.2 

42 

11.5 

27 

7.5 

18 

2.4 

43 

16.1 

5 

8.7 

9 

4.5 

44 

15.9 

66 

7.9 

9 

3.0 

45 

16.0 

~60 

6.3 

~37 

2.2 

46 

18.0 

7! 

6.5 

21 

3.7 

47 

17.7 

98 

6.2 

63 

3.3 

48 

9.7 

90 

-6.6 

11 

- 

49 

3.3 

31 

1.8 

8 

0.5 

50 

4.4 

93 

1.2 

4 

1.1 

51 

6.6 

94 

2.3 

30 

1.2 

*^max  =  maximum  in  the  test,  FJV  max  =  maximum  average  ice 
force,  Fit( j  =  standard  deviation  of  ice  force. 


45-51  the  reduction  of  averages  is  only  9-1 5%  but 
the  maximums  increase  13-21%. 

The  fourth  and  fifth  columns  of  Table  2  show 
the  maximums  of  the  averaged  ice  forces  and  their 
velocities  for  each  3.15-second  interval  of  data  col¬ 
lection.  This  time  period  is  long  enough  to  include 
several  ice  force  cycles  but  short  enough  to  show 
the  changes  with  increasing  velocity.  (The  graphic 
output  of  this  data  is  given  in  Appendix  A.)  The 
standard  deviations  and  the  maximum  values  during 
each  3.15-second  interval  are  given.  In  the  case  of 
the  maximum  of  the  averages  the  corresponding  ice 
velocity  is  always  low  and  matches  the  transition 
zone  from  ductile  to  brittle  loading  rate. 

An  apparent  discrepancy  with  the  theory  of  self- 
excited,  ice-induced  vibrations  is  that  the  pile  diam¬ 
eter  does  not  affect  the  velocity  of  the  occurrence 
of  maximums.  The  stress  rate  definition1  is 

8  ocv 
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where  ac  is  the  crushing  strength,  v  is  the  ice  velocity 
and  D  is  the  diameter  of  the  pile,  to  which  the  stress 
rate  is  inversely  proportional.  Because  the  stress 
rate  corresponding  to  the  transition  from  ductile  to 
brittle  should  be  the  same  in  the  same  ice  field,  the 
change  of  pile  diameter  from  64  to  1  78  mm  should 
have  increased  the  relevant  velocity  2.8  times.  There 
is  a  small  indication  of  this  trend  between  tests  37 
and  39,  but  elsewhere  the  results  are  independent. 


In  tests  50  and  51  the  initiation  of  buckling  prevents 
comparisons  to  test  49.  A  possible  explanation  for 
the  independence  of  the  pile  diameter  is  that  during 
pile  vibrations  the  relative  velocity  between  the  pile 
and  the  ice  is  lower  during  the  force  build-up  phase 
than  it  would  be  with  a  stiff  structure  with  ice  veloc¬ 
ity  alone.  Thus  the  same  stress  or  strain  rate  can  be 
achieved  in  a  wide  velocity  range,  or  inversely,  the 
same  stress  rate  can  be  achieved  with  the  same  veloc¬ 
ities  but  with  different  diameters. 

Another  factor  is  the  aspect  ratio  (pile  diameter 
at  waterline/ice  thickness)  effect.  With  larger  diam¬ 
eters  the  aspect  ratio  effect  is  smaller,  and  as  effective 
crushing  strength  should  be  used  in  eq  1,  it  tends  to 
reduce  the  diameter  effects.  If  this  is  taken  into  ac¬ 
count,  the  ratio  of  2.8  between  tests  37  and  39  re¬ 
duces  to  2.0.  The  questions,  "How  does  the  ice  in 
front  of  the  pile  feel  this  aspect  ratio  effect  and  how 
is  this  accounted  for  in  the  stress  rate  equation?"  are 
still  unanswered. 

Appendix  B  shows  that  the  standard  deviation  of 
the  ice  force  follows  quite  closely  the  behavior  of 
the  averages  with  increasing  velocities.  The  maximum 
standard  deviations  occur  in  the  area  of  transition 
from  ductile  to  brittle  ice  and  decrease  thereafter. 

The  ratio  of  the  average  to  the  standard  deviation  is 
usually  around  0.5,  but  it  varies  from  0.2  to  0.8.  In 
test  43  the  ratio  is  near  1 .0,  but  it  appears  that  the 
initial  zero  for  calculation  was  1.7  kN;  this  value 
should  be  added  to  both  the  average  and  the  maximum 


readings  in  that  test.  Thus,  the  correct  ratio  in  test 
43  is  around  0.47.  In  tests  34,  35  and  36  the  stan¬ 
dard  deviations  are  higher  than  average;  however, 
this  is  probably  a  result  of  errors  in  storing  data  on 
the  floppy  disk. 

If  a  short  sample  of  the  ice  force  history  is  avail¬ 
able,  the  maximum  ice  force  can  be  estimated  with 
the  following  equation 

'max^a v  +  3  •  Fs(d  .  (2) 

Equation  2  is  commonly  used  in  statistical  anal¬ 
ysis  in  predicting  maximums  for  normally  distributed 
samples.  A  further  study  is  planned  to  figure  out 
what  kind  of  statistical  distribution  the  measured 
ice  forces  follow.  The  coefficient  3  corresponds  to 
a  99.9%  confidence  limit.  Equation  2  predicts  the 
maximums  well  for  the  whole  velocity  range  in 
Appendix  B.  In  design,  therefore,  it  is  better  to 
estimate  the  maximum  ice  force  using  eq  2  than  to 
use  a  measured  maximum  value  directly. 

DYNAMIC  ASPECT  RATIO  EFFECT  AND 
CRUSHING  STRENGTH 

The  effective  ice  pressure  against  a  vertical  pile 
depends  on  the  pile  diameter  to  ice  thickness  ratio. 

A  great  many  reports  of  model  studies  and  in  a  few 
cases  of  full-scale  measurements  have  indicated  that 
the  ice  pressure  against  an  infinite  pile  has  to  be 


multiplied  by  a  factor  k  to  get  an  effective  pressure 
against  a  finite  width  pile:16 

k  =  +5H/D  (3) 

where  H  is  the  ice  thickness  and  D  is  the  pile  diam¬ 
eter  at  the  waterline.  Earlier  model  tests  have  been 
more  or  less  static  and  have  had  a  relatively  short 
stroke  of  the  test  piece  into  the  ice.  In  this  test 
series,  in  addition  to  varying  the  ice  thickness  and 
the  pile  diameter,  the  ice  velocity  was  varied  through¬ 
out  the  range  from  ductile  to  brittle  failure  mode,  and 
no  limitation  was  present  due  to  the  length  of  the 
ice  fields.  Thus,  the  dynamic  aspect  ratio  effect  dur¬ 
ing  continuous  crushing  could  be  measured. 

Table  3  gives  the  average  and  maximum  ice  forces 
and  the  effective  crushing  strengths  for  the  low  and 
high  velocity  ranges.  The  low  velocity  corresponds 
to  the  area  of  transition  from  ductile  to  brittle  ice, 
and  the  high  velocity  corresponds  to  the  brittle 
range,  with  the  possibility  that  the  pile  deflection 
rate  may  cause  ductile  loading  rates  during  force 
build-up.  For  better  comparisons  the  3.1 5-second 
samples  are  chosen  from  the  same  velocity  range  as 
much  as  possible.  Exceptions  are  due  to  the  small 
effect  of  longitudinal  cracks  and  in  tests  50  and  51 
to  buckling,  which  reduced  the  ice  loads  significantly 
in  the  low  velocity  range.  Computer  data  for  tests 
34-36  are  unreliable  for  velocities  above  20  mm/s. 

The  tests  for  each  ice  field  usually  look  1-2  hours, 
so  the  ice  in  the  last  tests  was  tempered  longer  and 
thus  was  weaker. 


Table  3.  Effective  crushing  strength  at  low  and  high  velocities. 


Low  velocity _  _ High  velocity 


Test 

number 

D 

(mm) 

H 

(mm) 

X* 

V 

(mm/s) 

(kN) 

^  max 
(kN) 

°av 

(kPa) 

°max 

(kPa) 

V 

(mm/s) 

(kN) 

F  max 
(kN) 

aav 

(kPa) 

°max 

(kPa) 

34 

102 

35 

1.65 

17 

8.2 

19.3 

1390 

3260 

90 

10.0 

1690 

35 

178 

38 

1.43 

15 

2.6 

7.1 

740 

90 

- 

13.2 

- 

1380 

36 

64 

38 

1  99 

17 

2.0 

4.5 

mm 

90 

- 

10.6 

- 

2230 

37 

64 

50 

2.22 

11 

4.2 

6.6 

El 

71 

2.6 

5.8 

370 

820 

38 

102 

51 

1.87 

14 

4.9 

8.5 

EH 

76 

3.3 

8.9 

340 

910 

39 

178 

52 

1.57 

11 

7.1 

14.5 

990 

73 

4.8 

1 1.6 

330 

800 

40 

178 

53 

1.58 

12 

9.4 

14.4 

- 

- 

- 

- 

- 

41 

178 

53 

1.58 

16 

9.0 

14.5 

75 

4.8 

12.2 

320 

820 

42 

102 

53 

14 

6.4 

10.7 

630 

1050 

76 

3.6 

10.0 

350 

980 

43 

102 

54 

1.91 

11 

6.7 

10.7 

640 

1020 

74 

3.4 

8.1 

320 

770 

44 

178 

55 

1.60 

9 

7.9 

12.8 

74 

4.9 

11.5 

740 

45 

102 

79 

2.21 

22 

5.7 

13.9 

780 

70 

5.4 

16.0 

900 

46 

178 

79 

20 

6.5 

13.3 

520 

70 

5.8 

710 

47 

178 

79 

1.80 

39 

5.3 

15.0 

240 

570 

63 

6.2 

16.3 

640 

48 

64 

82 

2.72 

24 

5.5 

7.9 

560 

90 

3.8 

9.7 

270 

690 

49 

64 

22 

1.65 

8 

1.8 

2.9 

780 

1260 

80 

2.6 

390 

1130 

50 

102 

21 

1.42 

4 

1.2 

3.2 

400 

1070 

80 

1.4 

3.5 

470 

1 180 

51 

178 

20 

1.24 

5 

1.3 

3.7 

860 

70 

2.1 

5.5 

490 

1280 

*  k  -  -v/ 1+5  H/D  ;  av  =  average  in  3.15-s  interval;  max  =  maximum  value;  o  =  T/kDH. 
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The  results  show  that  the  effective  ice  pressure 
can  be  predicted  quite  well  by  using  the  aspect  ratio 
effect  (eq  3).  The  calculated  crushing  strengths 
against  an  infinitely  wide  pile  given  in  Table  3  for 
different  diameters  are  relatively  close  to  their  aver¬ 
age  in  the  same  ice  field.  Even  the  results  from  ice 
fields  13  and  14  (tests  37-40  and  41-44)  are  almost 
identical  in  spite  of  the  dynamically  different  struc¬ 
tural  configurations.  Even  though  the  devia'  onsof 
crushing  stresses  in  each  ice  field  are  within  the  ex¬ 
pected  range,  there  is  a  trend  that  suggests  that  the 
aspect  ratio  effect  could  be  stronger  than  according 
to  eq  3. 

The  crushing  strength  decreased  with  increasing 
thickness  of  urea  ice.  The  explanation  is  that  the 
proportion  of  surface  layer  of  ice  with  a  vertical 
c-axis  is  much  larger  with  thin  ice  than  with  thick 
ice.  As  all  ice  fields  had  roughly  the  same  bending 
strength,  it  appears  that  the  surface  layer  does  not 
lose  its  crushing  strength  as  fast  as  it  lost  its  bending 
strength  during  tempering. 

In  these  tests  the  ratio  of  crushing  strength  to 
bending  strength  is  about  6  for  maximum  values 
and  about  3  for  average  values.  This  is  the  same 
ratio  as  for  full-scale  ice  strengths.  The  scale  reduc¬ 
tion  of  strength  due  to  urea-doping  and  tempering 
is  about  5  for  a  full-scale  flexural  strength  of  800 
kPa,  about  3  for  maximum  and  about  6  for  average 
crushing  strengths  if  a  value  of  2800  kPa  is  adopted 
for  the  full-scale  strength.  During  the  six-hour  tem¬ 
pering  the  flexural  strength  reduced  by  half,  the  rest 
of  the  reduction  being  due  to  urea.  No  tests  were 
conducted  to  determine  the  reduction  rate  of  the 
crushing  strength  during  tempering. 


MEASURED  ICE  FORCE  FREQUENCIES 

One  of  the  main  objectives  of  these  tests  was  to 
create  continuously  repeating  ice  force  fluctuations 
under  controlled  conditions.  In  most  tests  there  are 
good  examples  of  how  the  ice  force  repeats  itself 
with  similar  patterns  and  how  the  velocity  of  the  ice 
changes  both  the  shape  and  frequency  of  ice  force 
fluctuations.  In  test  40  the  ice  velocity  was  kept 
constant,  and  it  was  chosen  to  give  the  maximum 
amplitude.  A  repeating  force  and  displacement  pat¬ 
tern  developed  (Fig.  8a).  The  vibrating  pile  reached 
its  limit  cycles.  Only  the  randomness  in  ice  strength 
created  the  small  variations. 

The  frequency  of  the  ice  force  fluctuations  de¬ 
pended  on  the  ice  velocity  and  the  natural  frequencies 
of  the  structure.  The  effects  of  pile  diameter  and 
ice  thickness  were  also  present.  Table  4  includes  the 
dominant  frequencies  of  ice  forces  with  different 
velocities  in  those  tests  that  exhibited  distinct  ice 
force  fluctuations.  The  frequencies  were  measured 
directly  from  the  force-time  histories.  The  signal 
analyzer  was  not  used  to  solve  for  the  low  frequency 
components  because  both  long  sampling  periods  and 
several  averages  would  have  been  required,  which 
would  have  reduced  spectrum  peaks,  as  the  frequency 
was  continuously  changing  with  increasing  ice  veloc¬ 
ity.  Only  the  low  velocity  range  frequencies  are  given 
in  Table  4,  and  they  are  the  local  values  for  only  one 
or  two  loading  cycles  that  best  present  the  ice  force 
fluctuations  for  each  velocity  point. 

Plotter  graphs  (Appendix  C)  and  Table  4  indicate 
that  the  crushing  frequency  increases  at  the  same 
rate  as  the  ice  velocity  at  the  beginning  and  then 
starts  to  lag  at  velocities  over  20  mm/s.  Over  30-40 


Table  4.  Crushing  frequency  at  low  velocity  range.  During  each  test  three  samples  were  taken  as 
the  velocity  increased. 


Sample  1 _  _ Sample  2 _  _ Sample  3 


Fest 

number 

V* 

(mm/s) 

f 

(Hz) 

F 

(kN) 

f0 

(Hz) 

f/fq 

V 

(mm/s) 

f 

(Hz) 

F 

(kN) 

f0 

(Hz) 

f  t'o 

V 

(mm  Is) 

f 

(Hz) 

F 

(kN) 

f0 

(Hz) 

f/fp 

38 

10 

1.0 

6.7 

3.5 

0.29 

16 

2.1 

7.4 

5.0 

0.42 

27 

2.1 

8.1 

7.8 

0.27 

39 

10 

0.44 

14.6 

1.5 

0.29 

16 

1.3 

10.2 

3.7 

0.36 

27 

2.3 

10.3 

6.1 

0.38 

40 

12 

0.91 

14.4 

1.9 

0.47 

22 

2.4 

10.6 

5.0 

0.48 

- 

- 

- 

- 

41 

14 

0.96 

13.2 

2.5 

0.39 

19 

1.5 

14.1 

3.1 

0.48 

35 

2.1 

13.9 

5.9 

0.36 

42 

7 

0.67 

7.2 

2.3 

0.30 

18 

1.6 

10.1 

4.2 

0.48 

25 

2.0 

10.5 

5.6 

0.30 

43 

11 

0.82 

10.7 

2.4 

0.34 

19 

2.1 

7.8 

5.7 

0.39 

29 

2.2 

8.0 

8.5 

0.26 

44 

9 

0.48 

12.9 

1.6 

0.30 

17 

1.2 

12.9 

3.1 

0.37 

26 

2.3 

10.8 

5.6 

0.41 

Average 

10 

0.34 

18 

0.43 

28 

0.33 

v  -  ice  velocity,  f  -  measured  frequency, 

F  =  ice  force,  Ag  =  calculated  frequency  from  eq  5. 
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mm/s  the  crushing  frequency  jumps  intermittently 
to  8.5  Hz,  which  is  the  lowest  natural  frequency  of 
the  pile,  and  then  remains  at  8.5  Hz  as  the  velocity 
further  increases,  even  up  to  maximum  velocity 
(1 10  mm/s).  In  tests  43  and  44  the  velocity  range 
during  which  crushing  frequency  follows  ice  velocity 
was  below  20  mm/s;  above  that,  a  frequency  varying 
from  2.2  to  2.8  Hz  was  dominant  up  to  63  mm/s  in 
test  43  and  up  to  95  mm/s  in  test  44.  Thereafter 
the  crushing  frequency  increased  again  with  velocity. 
The  first  natural  frequency  of  the  pile  was  2.6  Hz 
and  that  of  the  carriage  was  2.5  Hz.  The  next  natural 
frequency  of  the  pile  in  tests  43  and  44  contributing 
significantly  to  the  deformation  of  the  pile  at  the 
waterline  was  37.5  Hz.  It  appears  that  the  natural 
frequencies  attract  crushing  frequencies.  This  proves 
that  the  response  of  the  structure  (in  addition  to  ice 
itself)  affected  the  resulting  ice  force. 

In  tests  38-44  the  ice  thickness  was  53±2  mm,  and 
the  stiffness  of  pile  in  the  direction  of  the  ice  force 
was  constant.  A  parameter  that  is  equally  as  impor¬ 
tant  for  crushing  frequency  as  ice  velocity  is  the 
magnitude  of  the  ice  force,  or  more  specifically,  the 
range  of  the  ice  force  during  repeating  loading  cycles. 
Table  4  verifies  that  crushing  frequency  is  inversely 
proportional  to  the  ice  force. 


CALCULATED  ICE  FORCE  FREQUENCIES 


There  are  two  reported  methods  for  predicting  the 
frequency  of  the  ice  force  during  crushing.  Neill16 
explained  that  ice  has  a  tendency  to  fail  into  floes 
of  a  certain  size,  0.2-0.5  times  the  ice  thickness. 

This  yields  the  equation 


f=  2 


(4) 


Another  method12  starts  from  the  length  of  ice  that 
is  crushed  during  each  cycle.  If  it  is  equal  to  the 
elastic  displacement  of  the  structure  with  the  max¬ 
imum  ice  force  just  before  crushing  starts,  then 


f  = 


ki < 

T 


(5) 


where  k  is  the  stiffness  coefficient  of  the  structure 
toward  the  direction  of  ice  action  at  the  ice  action 
point.  If  the  ice  forces  do  not  reduce  to  zero  during 
the  crushing  phase,  the  range  of  the  ice  force  should 
be  used  instead  of  the  maximum.  In  dynamic  con¬ 
ditions  the  magnification  factor  should  be  applied 
to  the  ice  force  to  account  for  the  increasing  crushing 
lengths. 

If  eq  4  is  used  to  predict  the  low  crushing  fre¬ 
quencies  in  Table  4,  it  succeeds  in  all  but  three  of 
the  intermediate  velocity  cases.  However,  it  falls 


short  of  predicting  the  constant  frequency  of  8.5 
Hz  (all  the  tests  except  43  and  44).  In  tests  43  and 
44  the  frequency  of  2.5  Hz  can  be  predicted,  but 
only  with  an  unrealistic  assumption  that  the  coeffic¬ 
ient  0.2-0.5  in  the  denominator  should  be  velocity 
dependent.  In  the  tests  the  crushed  or  failed  ice 
pieces  were  never  larger  with  higher  velocities.  Also, 
in  tests  43  and  44  eq  4  fails  to  predict  crushing  fre¬ 
quencies  over  15  Hz. 

The  columns  for  fQ  in  Table  4  present  the  crushing 
frequencies  calculated  using  eq  5.  It  appears  that  the 
prediction  is  always  too  high.  It  should  be  romewhat 
higher  because  eq  5  takes  into  account  only  the  time 
for  load  build-up,  not  the  crushing  phase.  However, 
from  plotter  graphs  it  can  be  seen  that  the  crushing 
phase  is  very  short,  and  thus  its  contribution  to  crush¬ 
ing  frequency  is  negligible  at  low  velocities.  Another 
factor  is  that  the  inertia  of  the  structure  may  make 
the  crushed  distance  longer  than  the  static  displace¬ 
ment  Fjk  in  eq  5. 

The  ratio  of  measured  to  calculated  frequencies 
averaged  from  0.33  to  0.43  with  relatively  small  de¬ 
viations.  As  both  velocity  and  ice  force  in  eq  5  arc 
measured  quantities,  a  reduction  coefficient  should 
be  applied  to  the  pile  stiffness  6.  As  ice  is  pushed  by 
the  test  carrier  bridge,  the  stiffness  of  the  bridge, 
frames,  ice  boom  and  ice  also  have  to  be  taken  into 
account.  These  stiffnesses  arc  17.2,  155,  1 9.5  and 
1 1 2  kN/mm,  respectively.  In  series  with  the  pile 
stiffness  of  2.33  kN/mm,  the  overall  stiffness  is  1 .80 
kN/mm,  yielding  a  reduction  coefficient  of  0.77.  In 
addition,  gaps  in  the  ice  boom  joints  and  frame  con¬ 
nections  decrease  the  stiffness  further. 

The  effective  stiffness  of  the  pile  against  the  ice 
force  can  be  calculated  from  plotter  graphs  using  the 
definition  of  stiffness: 


F  _  AF  _  A  F 
x  Ax  v-At 


(6) 


During  the  load  build-up  phase  a  sample  taken  over 
a  time  of  At  is  chosen  that  yields  a  displacement 
vAt  and  an  increase  in  ice  force  AF.  Table  5  pro¬ 
vides  data  for  calculating  effective  stiffnesses  in  some 
tests. 

Straight  load  build-up  in  plotter  graphs  exist  only 
at  low  velocities.  The  loading  rate  is  then  so  low  that 
viscoelastic  yielding  increases  the  deformations  in  ice 
and  thus  decreases  the  calculated  stiffness  values,  be¬ 
cause  velocity  v  in  Table  5  is  ice  velocity  and  not  the 
relative  velocity  between  the  ice  and  the  structure. 
Another  factor  is  the  contact  coefficient.  The  less 
area  of  pile  that  touches  the  ice,  the  higher  the  local 
ice  pressures  and  displacements  and  the  lower  the 
effective  stiffness. 

The  frequency  ratios  in  Table  4  coincide  with  the 
stiffness  ratios  in  Table  5,  recognizing  the  viscoelastic 
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Table  5.  Calculation  of  effective  stiffnesses. 


Test 

number 

At 

(S) 

V 

(mm/s) 

AF 

(kN) 

weff 

(N/mm) 

^effl^olle 

39 

3.0 

4.6 

8.0 

580 

0.25 

39 

2.0 

9.0 

15.0 

840 

0.36 

40 

3.3 

3.9 

10.3 

800 

0.34 

41 

3.0 

9.7 

11.7 

400 

0.17 

42 

3.0 

4.2 

6.5 

520 

0.22 

42 

2.0 

8.8 

7.2 

410 

0.18 

43 

3.0 

4.2 

11.1 

880 

0.38 

44 

2.0 

4.5 

5.9 

660 

0.28 

44 

1.0 

8.9 

4.5 

500 

0.21 

44 

2.0 

4.6 

6.6 

720 

0.31 

Average 

630 

0.27 

effects  in  the  latter.  The  measurements  of  frequen¬ 
cies  correlate  well  with  the  values  predicted  by  eq 
5;  the  correct  frequencies  can  be  predicted  by  using 
the  effective  stiffness  instead  of  the  pile  stiffness 
alone.  In  addition  the  frequency  ratios  according 
to  eq  5  are  lowest  at  low  loading  rates,  where  visco¬ 
elastic  effects  are  pronounced;  they  reach  a  maximum 
as  load  build-up  time  shortens  and  then  decrease  at 
increasing  velocities  as  the  time  for  deflection  spring- 
back  becomes  more  significant.  This  measured  de¬ 
pendence  was  predicted  earlier  using  computer  sim¬ 
ulations.15 


a.  Tests  37-42. 


Equation  5  also  neatly  explains  why  the  natural 
frequencies  of  the  structure  can  persist  in  the  ice 
force  for  a  wide  velocity  range.  The  range  of  ice 
force  and  the  length  of  the  crushed  ice  are  increasing 
with  velocity,  while  the  frequency  can  remain  con¬ 
stant.  For  example,  in  test  41  the  first  natural  fre¬ 
quency  of  the  pile  (8.5  Hz)  appears  at  40  mm/s  veloc¬ 
ity  and  persists  up  to  the  maximum  velocity  of  85 
mm/s.  The  range  of  the  ice  force  without  dynamic 
correction  at  40  mm/s  is  4  kN  and  at  85  mm/s,  8  kN. 
If  the  effective  stiffness  is  850  N/mm,  the  crushing 
length  increases  from  4.7  to  9.4  mm. 

Peyton14  measured  ice  force  frequencies  at  Cook 
Inlet  with  a  full-scale  structure  that  had  the  first 
natural  frequency  around  1.0  Hz.  Because  measure¬ 
ments  under  varying  conditions  yielded  a  crushing 
frequency  of  about  1  Hz,  he  concluded  that  it  is  a 
physical  property  of  ice  to  fail  at  this  frequency  dur¬ 
ing  crushing.  In  the  present  tests,  8.5-Hz  vibrations 
also  prevailed  at  varying  ice  thicknesses  and  velocities. 
However,  when  the  structural  configuration  was 
changed  and  there  was  no  dynamically  unstable  nat¬ 
ural  mode  at  8.5  Hz,  ice-induced  vibrations  no  longer 
appeared  at  this  frequency.  This  proves  that  the 
crushing  frequency  is  not  a  property  of  ice  alone 
but  is  also  the  result  of  the  structural  properties  in 
a  dynamic  ice-structure  interaction. 

A  signal  analyzer  was  used  to  calculate  the  8.5-Hz 
vibration  amplitudes  (Fig.  10).  In  tests  37-42  the 


b.  Tests  46-51 . 


Figure  10.  Ice  force  range  at  8.5  Hz  us  velocity  (ice  force  without  dynamic  correction). 
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pile  diameter  had  no  apparent  effect  on  the  vibra¬ 
tion  amplitudes.  Theoretically  when  stress  rate  ef¬ 
fect  is  taken  into  account,  there  is  no  direct  depend¬ 
ence  on  diameter  either. 

Substituting  the  velocity  fr- — *  eq  1  and  the  ice 
force  components  into  eq  5  results  in 

.  7 tDo 

*  0 

Skoc  *  k  o 

f  =  - =  1  .  — 1 —  —  (7) 

kocHD  8  K2acH  ac 

which  shows  that  the  diameter  has  its  effect  only 
through  the  square  of  the  aspect  ratio  coefficient  k. 

The  aspect  ratio  effect  k  is  taken  into  account  also 
in  the  stress  rate  (eq  1 ),  because  oc  in  eq  1  is  the 
maximum  elastic  stress  in  the  ice  in  front  of  the 
pile.  In  tests  4 5-51  the  ice  thickness  was  changed, 
and  Fig.  Mb  shows  that  there  is  not  much  change 
in  amplitudes  compared  to  tests  37-42.  The  reason 
is  that  the  crushing  strength  is  lower  with  thicker 
ice  and  higher  with  thin  ice  than  in  tests  37-42 
(Table  3).  This,  together  with  the  aspect  ratio  ef¬ 
fects,  results  in  small  changes  in  crushing  frequencies, 
or  inversely,  in  crushing  amplitudes,  because  the 
frequency  is  constant  because  the  natural  frequency 
of  the  structure  dominates. 

A  significant  change  in  crushing  amplitudes  oc¬ 
curred  in  tests  49  and  50.  With  thin  ice  the  ice  force 
fluctuations  were  small,  and  when  the  velocity  was 
high  enough,  eq  5  predicted  higher  frequencies  than 
the  natural  frequency  of  8.5  Hz.  With  increasing 
velocity  the  amplitudes  first  increased  as  long  as  the 
undisturbed  frequency  remained  below  8.5  Hz  and 
reduced  thereafter.  The  undisturbed  frequency  is 
the  crushing  frequency  that  would  prevail  if  the 
natural  frequency  of  the  pile  did  not  dominate  the 
response;  compare  the  low  velocity  frequencies  in 
tests  43  and  44  (Table  4).  Thus,  as  frequency  is 
constant,  a  rcsonancc-likc  condition  appears  in  the 
amplitude  response.  In  test  50  the  diameter  was 
wider,  requiring  a  higher  velocity  for  resonance.  For 
test  5  I  the  resonant  velocity  was  beyond  the  lest 
velocity  range.  With  the  thick  ice  in  tests  37-42  and 
45-48  the  velocity  never  reached  resonance  with 
8.5  Hz. 

Aj  the  velocity  increased  over  84  mm/s  in  test 
49  and  over  92  mm/s  in  test  50,  the  second  natural 
frequency  (60  Hz,  bru  .not  AC  noise;  compare  Fig. 

4)  started  to  appear  in  the  ice  force.  For  short  per¬ 
iods  the  first  mode  disappeared  and  the  response 
was  almost  entirely  second  mode.  Again  eq  5  gives 
the  correct  frequency,  with  the  measured  ice  force 
range  of  I  kN  and  effective  pile  stiffness  of  850  N/mm. 


In  tests  43  and  44  the  first  two  natural  frequen¬ 
cies  of  the  pile  were  2.6  and  37.5  Hz.  According  to 
stability  theories  for  ice-induced  vibrations,  the  2.6- 
Hz  mode  should  be  stable  and  should  not  be  excited 
by  the  ice  crushing.  The  averaged  power  spectrums 
indicate  that  this  is  the  case  in  test  43.  There  is  some 
excitation  in  the  first  superstructure  mode  up  to  100 
mm/s,  but  resonance  never  occurs. 

In  test  44  the  only  difference  was  the  pile  diam¬ 
eter,  which  only  increases  the  ice  forces  when  the 
dynamic  behavior  of  the  pile  remains  unaltered.  In 
this  test,  however,  continuous  high-level  vibrations 
occurred  in  a  frequency  range  of  2. 2-2.8  Hz,  while 
the  ice  velocity  was  between  35  and  75  mm/s.  The 
most  persistent  vibrations  occurred  at  2.5  Hz,  and 
the  maximum  amplitudes  were  recorded  at  a  velocity 
of  60  mm/s,  which  can  be  regarded  as  a  resonant 
point.  The  analysis  of  ice  force  and  acceleration 
waveforms  revealed  that  it  was  not  a  resonance  with 
the  lest  pile,  since  the  phase  shift  with  the  ice  force 
and  pile  acceleration  is  1 42°;  more  likely  it  is  a  res¬ 
onance  with  the  carriage,  for  which  the  phase  shift 
is  112°  (Fig.  11).  This  condition  prevailed  in  the 
whole  velocity  range  from  35  to  75  mm/s  for  which 
2.5-Hz  vibrations  were  continuous.  In  tests  for  the 
carriage  alone  the  natural  mode  with  a  frequency  of 
2.5  Hz  was  measured.  It  appears  that  the  flexible 
carriage  and  ice  field  is  a  mechanical  system  that  is 
excited  by  the  ice  crushing.  This  system  is  dynamic¬ 
ally  unstable  and  is  the  reason  for  the  resonance 
around  2.5  Hz.  As  the  first  natural  frequencies  of 
the  test  pile  and  carriage  were  almost  identical,  the 
vibrations  of  the  test  pile  were  also  strongly  ampli¬ 
fied. 

The  dependence  of  crushing  frequency  and  ice 
force  amplitudes  on  the  whole  velocity  range  in  these 
tests  can  be  explained  using  eq  5.  A  similar  depend¬ 
ence  can  also  be  predicted  analytically  it  the  dynamic 
response  of  the  pile  is  lied  together  with  the  depend¬ 
ence  of  the  ice  crushing  on  the  loading  rale.  Thus, 
it  can  be  concluded  that  eq  5  is  valuable  in  predicting 
crushing  frequency  ranges  and  possibilities  lor  res¬ 
onance.  The  natural  frequencies  ol  the  structure 
persist  in  a  wide  velocity  range,  while  force  and  dis¬ 
placement  amplitudes  change  instead  of  frequency. 
Otherwise  the  frequencies  that  eq  5  predicts  arc 
always  a  bit  too  high.  This  is  due  to  viscoelastic 
yielding  at  low  velocities,  deflection  spring-back 
times  at  high  velocities,  and  reduced  ellcctive  pile- 
carriage  stiffness.  In  the  field  the  latter  results  only 
Irom  the  elastic  deformations  ol  ice  and  is  due  to 
uneven  contact  between  the  ice  and  the  pile. 
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a.  Test  pile  top  acceleration  phase  vs 
carriage  acceleration,  and  carriage  accel¬ 
eration  amplitude  spectrum. 
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b.  Carriage  acceleration  phase  its  ice 
force,  and  carriage  acceleration  ampli¬ 
tude  spectrum. 


Ice  Force 
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c.  Phase  vector  diagram  for  ice  force,  pile  and  carriage 
accelerations. 


Figure  1 1.  Ice  force,  pile  and  carriage  acceleration  phases  ( test  44,  v  -  30-60  mm/s). 
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Table  6.  Accelerations,  velocities  and  displacements. 


Acceleration  Velocity f  Displacement/ 


Test 

''ice 

(cm/s) 

(9)* 

(cm/f) 

(cm) 

Note 

parallel 

perpendicular  parallel  perpendicular 

parallel 

perpendicular 

37 

2.0 

1.2 

0.4 

22 

7 

0.41 

0.14 

Decaying 

10.1 

3.2 

2.2 

59 

- 

1.11 

- 

Continuous 

2.0 

2.6 

0.4 

48 

6 

0.11 

0.14 

Decaying 

38 

8.9 

2.4 

1.0 

45 

18 

0.83 

0.34 

Continuous 

10.9 

2.6 

1.0 

48 

18 

0.11 

0.34 

Continuous 

39 

1.5 

2.0 

0.5 

37 

9 

0.70 

0.17 

Decaying 

8.7 

2.0 

- 

37 

- 

0.70 

- 

Continuous 

40 

1.2 

4.0 

0.7 

74 

12 

1.39 

0.22 

Decaying 

1.2 

3.0 

2.2 

56 

41 

1.04 

0.76 

Continuous 

41 

1.8 

3.6 

0.5 

67 

9 

1.25 

0.17 

Decaying 

8.7 

2.4 

2.1 

45 

39 

0.83 

0.73 

Decaying 

42 

1.7 

7.7 

0.7 

- 

12 

- 

0.22 

Decaying 

6.4 

1.6 

0.4 

30 

8 

0.56 

0.14 

Continuous 

1.2 

8.2 

2.1 

- 

16 

_ 

0.12 

Decaying 

43** 

3.3 

8.5 

2.6 

- 

19 

- 

0.15 

Decaying 

9.5 

8.9 

1.8 

66 

14 

0.51 

0.10 

Continuous 

1.0 

8.3 

3.0 

_ 

» 

_ 

44** 

2.7 

7.1 

1.6 

- 

12 

- 

0.09 

Decaying 

6.7 

9.3 

3.1 

- 

- 

- 

- 

*1 

6.5 

0.9 

* 

_ 

_ 

45 

*6 

8.5 

0.7 

- 

- 

- 

- 

=*9 

<2.8 

0.8 

<51 

- 

<0.96 

- 

2.1 

8.5 

1.3 

_ 

*1 

_ 

46 

7.2 

8.9 

2.6 

<81 

- 

<1.5 

- 

7.2 

1.5 

- 

28 

- 

0.52 

- 

Continuous 

1.7 

2.0 

2.0 

_ 

_ 

_ 

47 

5.0 

8.5 

1.6 

- 

- 

- 

- 

9.7 

2.0 

*2.0 

37 

- 

0.70 

- 

^Continuous 

9.0 

2.4 

_ 

44 

0.83 

_ 

^Continuous 

48 

9.8 

2.4 

- 

44 

- 

0.83 

- 

Continuous 

9.8 

2.6 

1.3 

48 

- 

0.90 

- 

Continuous 

2.5 

0.6 

0.25 

1  1 

4.5 

0.20 

0.08 

Continuous 

49 

6.9 

1.7 

- 

31 

- 

0.58 

- 

Continuous 

8.8 

2.8 

- 

.  7.1 

- 

0.02 

- 

Continuous 

2.2 

0.6 

- 

11 

_ 

0.20 

_ 

Continuous 

50 

7.8 

1.4 

- 

26 

- 

0.48 

- 

Continuous 

9.2 

2.4 

- 

6.1 

- 

0.02 

- 

Continuous 

51 

2.0 

4.2 

0.7 

_ 

_ 

_ 

9.3 

2.4 

1.1 

44 

- 

0.8  1 

- 

Continuous 

•Acceleration  in  units  of  gravity  of  earth  (9.8  I  m/s^(. 

fVelocity  and  displacement  integrations  were  carried  out  only  when  a  distinct  mode  was  present. 

••In  tests  43  and  44  longitudinal  accelerations,  velocities  and  dicplacements  were  measured  at  the  top  ot 
lower  structure  and  at  the  bottom  of  superstructure  instead  of  parallel  and  perpendicular  to  the  ice  torce. 


ACCELERATIONS,  VELOCITIES  AND 
DISPLACEMENTS 


The  acceleration  transducers  at  the  superstructure 
were  used  to  measure  longitudinal  and  cross-axis 
vibrations  in  all  but  tests  43  and  44,  where  they 
measured  the  longitudinal  accelerations  just  above 
and  below  the  vibration  isolation  part.  This  meas¬ 
ured  the  effectiveness  of  the  vibration  isolation. 
Table  6  contains  the  data  on  maximum  accelerations 
during  each  test  at  different  velocities.  If  the  accel¬ 
eration  waveform  is  distinguishable  as  a  single  mode, 
then  velocities  and  displacements  were  integrated. 
For  cases  of  decaying  vibrations  the  initial  maximum 
values  are  given. 

During  the  low  velocity  load  build-up  phase  the 
acceleration  was  zero.  During  ice  crushing,  high  fre¬ 
quencies  from  81  to  1 1 8  Hz  were  present.  These 
are  higher  than  the  two  lowest  natural  frequencies 
of  the  structure.  H igh  frequency  components  die 
out  fast,  and  simultaneously  the  first  or  second 
natural  modes  become  dominant.  The  crushing  is 
like  an  impact  loading  that  excites  all  the  natural 
frequencies  of  the  structure,  but  only  the  lowest 
exhibit  significant  displacements.  After  crushing, 
the  first  mode  vibrations  die  out  exponentially. 

This  first  mode  frequency  is  significantly  higher 
than  the  free  natural  mode  of  the  pile.  As  the  ice 
edge  sticks  to  the  pile  after  crushing,  it  induces  a 
supporting  effect,  which  increased  the  natural  fre¬ 
quency  in  these  tests  from  8.5  Hat  to  1 2-1  5  Hz. 

With  higher  velocities,  when  the  first  natural 
mode  of  the  structure  is  able  to  dominate  the  crush¬ 
ing  frequencies,  the  main  waveform  of  acceleration 
is  a  8.5-Hz  sinusoidal  curve  with  a  second  mode  at 
60  Hz  superimposed  on  it.  During  the  crushing 
phase  even  higher  components  arc  present.  These 
induce  the  highest  momentary  accelerations.  The 
maximum  recorded  values  are  over  9  g.  In  reality 
they  have  been  higher,  since  1 )  the  maximum  range 
of  the  accelerometers  was  exceeded,  which  caused 
saturation,  2)  these  peaks  arc  also  due  to  frequencies 
above  100  Hz  and  the  accelerometers  had  natural 
frequencies  around  100  Hz  with  60%  critical  damp¬ 
ing,  which  has  an  ef  fect  of  a  low  band  pass  filter, 

3)  the  data  were  recorded  on  a  tape  recorder  with 
a  0-1  56  Hz  hand,  and  4)  the  plotter  attenuates 
amplitude  peaks  over  50  Hz. 

Transverse  vibration  peaks  may  be  as  high  as  lon¬ 
gitudinal  ones  but  generally  they  were  less  than  one 
thiul  as  high.  Randomness  is  the  main  cause  of 
transverse  excitations  while  ice  is  failing  first  on 
one  side  at  the  pile  and  then  the  other.  If  contin¬ 
uous  longitudinal  vibrations  persist,  the  transverse 
vibration  remains  low. 

I  he  integrated  velocities  at  the  top  of  the  pile 


verify  that  the  deformation  rate  of  the  pile  can 
easily  be  much  higher  than  the  ice  velocity.  Consid¬ 
ering  the  shape  of  the  first  mode  vibrations,  ampli¬ 
tudes  at  the  top  should  be  multiplied  by  0.238  to 
get  those  at  the  ice  action  point.  Thus,  a  ductile 
loading  rate  is  possible  over  the  whole  velocity  range. 
The  maximum  velocities  are  800  mm/s  at  the  top  and 
1 90  mm/s  at  the  waterline.  The  maximum  integrated 
top  displacements  are  14  mm,  which,  however,  have 
been  dampened. 

In  tests  43  and  44  the  vibration  isolation  part 
dampened  the  vibrations  in  the  superstructure  to 
one  third  or  one  fourth  of  those  at  the  lower  struc¬ 
tures.  During  the  carriage  resonance  part  of  test  44 
(2.6  Hz)  the  acceleration  amplitude  in  the  super¬ 
structure  was  0.82  g.  This  resulted  in  a  velocity  of 
493  mm/s  and  a  displacement  of  30  mm. 

In  tests  41-43  one  longitudinal  accelerometer  was 
on  the  ice;  in  test  44  it  was  on  the  carriage.  During 


T ime  (ms) 


a.  Ice  in-plane  longitudinal  acceleration. 


OO  80  1 20  160 

Time  ( ms ) 


b.  Carriage  longitudinal  acceleration  at  84  Hz. 
Figure  12.  Acceleration  during  ice  crushing. 
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Table  7.  Damping  during  crushing. 


the  crushing  phase,  high-frequency  vibrations  over 
80  Hz  could  be  measured  on  both  the  ice  and  the 
carriage.  The  maximum  recorded  accelerations  were 
0.68  and  0.52  g,  respectively.  As  previously  men¬ 
tioned,  these  peaks  may  have  been  higher  but  they 
were  not  recorded  due  to  the  limits  of  the  instru¬ 
mentation.  Continuous  vibration  in  the  ice  with  2.6  Hz 
caused  an  acceleration  of  0.04  g,  a  velocity  of  24 
mm/s  and  displacement  amplitudes  of  1 .4  mm.  In 
the  carriage  the  corresponding  figures  were  0.08  g, 

48  mm/s  and  2.9  mm. 

Natural  modes  for  80-1 18  Hz  vibrations  enoun- 
tered  both  in  the  ice  and  at  the  carriage  cannot  be 
accounted  for.  Not  enough  accelerometers  were 
installed  for  identifying  these  modes.  Either  the 
natural  modes  of  the  carriage,  the  ice  field  in-plane 
vibrations,  or  the  combination  of  these  two  were 
excited  during  the  crushing  phase. 

Samples  of  high-frequency  vibrations  in  the  ice 
in  test  43  and  at  the  carriage  in  test  44  are  presented 
in  Fig.  12.  In  both  cases  vibration  with  a  frequency 
around  85  Hz  has  started  almost  at  the  same  time  as 
crushing;  thereafter  the  amplitudes  remained  at  a 
constant  level  for  as  long  as  crushing  lasted.  This 
might  be  because  the  mode  in  question  is  unstable 
for  ice-induced  vibrations.  Hence,  limit  cycles  de¬ 
velop  fast,  and  thereafter,  only  random  variations 
in  vibration  amplitudes  occur. 

The  frequency  range  from  80  to  1 18  Hz  looks 
high,  but  when  these  scale-model  tests  are  converted 
to  full  scale,  the  resulting  frequency  range  is  around 
10  Hz.  This  is  below  the  1 5-Hz  vibrations  measured 
from  a  bridge  pier  at  Pembina  River.19  With  a  1 :8 
scale  factor  these  model  tests  and  the  measurements 
at  Pembina  River  are  comparable.  The  ice  velocity 
at  Pembina  River  was  higher  at  the  beginning  of 
crushing  but  then  decreased  to  zero  as  ice  floes  in 
the  river  were  stopped  by  the  pier.  The  simply  sup¬ 
ported  pier  had  a  lowest  natural  frequency  of  15 
Hz  and  had  low  internal  damping,  which  evidently 
made  it  unstable  for  ice-induced  self-excited  vibra¬ 
tions. 


DAMPING 

The  damping  of  the  test  pile  was  measured  both 
in  and  out  of  the  water,  for  the  first  lowest  natural 
modes  there  was  no  sigml  leant  change  in  dantpin; 
or  in  frequencies.  At  8.5  Hz,  damping  was  0.8/ 
from  the  critical  damping,  and  at  60.0  Hz  it  was 
2.7%.  For  tests  43  and  44  damping  at  2.6  Hz  was 
3.2%,  and  at  37.5  Hz  it  was  7.9%.  Damping  at  2.6 
and  8.5  Hz  was  measured  from  the  decay  rate  of 
free  oscillations;  at  other  frequencies  it  was  measured 


Test 

number 

A/* 

A„ 

Number 
of  cycles 

r 

(96) 

34 

455 

100 

7 

3.4 

34 

365 

84 

7 

3.3 

35 

158 

36 

3 

7.9 

35 

82 

21 

3 

7.2 

36 

90 

30 

3 

5.8 

39 

150 

38 

4 

5.5 

40 

245 

43 

3 

9.2 

40 

210 

43 

6 

4.2 

40 

210 

42 

6 

4.3 

40 

190 

28 

6 

5.1 

41 

219 

58 

5 

4.2 

45 

215 

115 

2 

5.0 

45 

290 

155 

2 

5.0 

46 

170 

60 

3 

55 

46 

316 

72 

4 

5.9 

47 

263 

152 

3 

2.9 

Average 

5.3 

*A  ,  =  initial  amplitude,  An  =  amplitude 
after  n  cycles,  jr  =  damping  factor  (A  = 

2 7r  ?  =  logarithmic  decrement). 


from  autospectrum  half-power  points  in  conjunc¬ 
tion  with  transfer  function  measurements  using 
the  signal  analyzer.  The  windowing  effect  and  weight¬ 
ing  functions,  together  with  averaging  in  the  signal 
analyzer,  broaden  the  peaks,  resulting  in  considerably 
higher  damping  ratios  than  are  measured  directly  from 
the  decaying  signal.  Thus  the  signal  analyzer  damping 
ratios  should  be  used  only  for  comparisons. 

The  amount  of  damping  during  the  tests  can  be 
determined  cither  from  the  width  oi  the  peaks  of  the 
signal  autospcctrums  or  from  the  decay  rates.  Table 
7  presents  the  measured  dampings  during  the  decaying 
vibration  after  crushing  in  low  velocity  range.  The 
results  show  that  there  was  a  significant  increase  in 
damping  (from  0.8  to  5.3%)  due  to  the  ice  presence. 
As  the  main  crushing  had  already  occurred,  the  damp¬ 
ing  mechanism  must  have  been  due  to  friction  be¬ 
tween  broken  ice  pieces  and  the  pile  and  to  local 
crushing  of  peaks  in  the  uneven  ice  edge  while  a  new 
contact  developed. 

There  was  randomness  in  the  amount  ol  damping. 

In  some  cases  there  was  no  significant  vibration  after 
uushing,  which  means  that  all  the  elastic  energy  of 
the  pile  had  been  dissipated  during  crushing.  Similar 
phenomena  have  been  observed  in  the  field  with 
flexible  structures.  It  is  evident  that  the  stiller  the 
structure,  the  less  chance  there  is  for  damping  be¬ 
cause  of  ice  being  ground  into  small  pieces  between 
the  pile  and  the  ice  edge.  To  apply  this  measured 
damping  increase  to  full-scale  structures,  one  needs 
to  observe  the  displacement  amplitudes  carefully. 


ICE-INDUCLD  NtGATIVt  DAMPING 


With  higher  velocities  during  resonance  at  2.6  Hz 
in  test  44  the  damping  measured  from  the  ice  force 
autospectrum  peak  width  indicates  a  decrease  from 
19  to  7.8%.  This  verifies  the  concept  of  ice-induced 
negative  damping,  which  follows  theoretically  from 
the  negative  slope  in  the  ice  crushing  strength  versus 
loading  rate  curve. 

The  measurement  of  negative  slopes  from  com¬ 
puter  plots  (Appendix  A)  did  not  succeed  for  all 
tests.  Table  8  contains  the  data  from  those  tests 
where  negative  slopes  for  average  ice  forces  appeared. 
The  best  examples  are  tests  43  and  49,  followed  by 
tests  37  and  39.  In  Table  7,  v  is  the  ice  velocity  and 
not  the  relative  velocity  between  the  ice  and  pile. 
However,  at  the  end  of  the  load  build-up  phase  as 
the  pile  movement  stops  and  crushing  starts,  the 
relative  velocity  is  the  same  as  the  ice  velocity.  The 
damping  factor  f  is  for  a  model  with  a  single  degree 
of  freedom: 


dF 

du 


CO 

2k 


100 


(8) 


where  co  is  pile  vibration  angular  frequency  and  k 
is  pile  stiffness.  Equation  8  results  directly  from 
the  definition  of  damping  force. 

In  most  tests  (and  evidently  in  all  tests  at  some 
velocity  intervals)  there  is  a  negative  slope  so  steep 
that  it  induces  aperiodic  divergence  for  ice-induced 
vibrations.  This  occurs  as  crushing  starts,  after  which 
the  pile  accelerates  rapidly  to  so  high  a  velocity  that 
the  relative  velocity  between  the  ice  and  the  pile  cor¬ 
responds  to  a  velocity  in  the  brittle  region,  with  an 
average  zero  slope  in  the  ice  strength  curve.  With 
higher  modes  the  negative  damping  coefficient  is 
still  larger,  which  explains  why  the  ice  action  point 
of  the  structure  can  accelerate  so  rapidly  and  make 
the  crushing  phase  short. 


Table  8.  Negative  damping 
coefficients. 


rest 

number 

dF/3v 

(kNsIcm) 

r 

<%) 

*r 

(%) 

37 

-1.3 

-150 

9.8 

38 

-0.7 

-  70 

10.5 

39 

-1.7 

-200 

10.6 

41 

-1.4 

-160 

10.6 

42 

-1.4 

-160 

10.4 

43 

-3.6 

-130 

9.2 

44 

-0.7 

-  30 

7.8 

49 

-0.45 

-  50 

10.3 

I  he  last  column  in  I  able  7  gives  the  damping 
ratios  at  8.5  Hz  (or  2.6  Hz  in  tests  43  and  44)  during 
continuous  crushing.  These  ratios  are  measured  from 
the  widths  of  ice  force  autospectra.  The  ratio  fr  is 
the  equivalent  viscous  damping,  which  averages  the 
real  nonlinear  damping  during  the  whole  vibration 
cycle.  Thus,  even  though  ice-induced  negative  damp¬ 
ing  is  supercritical,  it  appears  oniy  for  a  part  of  the 
vibration  cycle;  the  result  is  that  damping  ratios  re¬ 
mained  practically  unchanged  at  8.5  Hz  and  decreased 
from  19.0  to  7.8%  at  2.6  Hz.  The  overall  effect  of  ice- 
induced  negative  damping  is  more  than  these  figures 
indicate  because  the  ice  also  induces  positive  damping. 
This  contribution  should  be  added  if  the  total  nega¬ 
tive  damping  effect  is  desired. 


LIMIT  CYCLES 

Measurement  data  plots  show  that  the  ice  force 
repeats  itself  continuously  in  consistent  patterns. 

The  vibrations  of  the  pile  have  reached  a  steady  state 
limit  cycle.  This  means  that  autonomous  self-excited 
vibrations  have  developed.  The  reason  for  the  rise  of 
vibrations  is  the  high  negative  ice-induced  damping, 
which  causes  dynamic  instability  in  the  small  (the 
growth  of  infinitesimal  disturbances  into  distinguish¬ 
able  increasing  dynamic  response).  However,  stability 
in  the  large  exists  because  vibrations  and  ice  forces 
fluctuations  stabilize  at  a  steady  level. 

The  limit  cycles  can  be  predicted  analytically  by 
starting  with  the  averaged  ice  crushing  strength  curves 
and  taking  into  account  the  dynamic  response  of  the 
structure,  yielding  similar  response  and  ice  force  pat¬ 
terns  to  those  that  were  measured.  Figure  13  is  the 
numerically  integrated  ice  force  history  for  test  40. 
Even  though  there  are  differences,  the  force  pattern 
and  frequency  correspond  to  the  measured  response 
(Fig.  8a).  Analytically  predicted  stable  and  unstable 
modes  were  verified  by  the  measurements.  The  8.5- 
Hz  mode  of  the  test  pile  was  very  unstable,  and  the 
measurements  exhibited  continuous  vibrations.  The 


*f  =  momentary  damping  factor, 
Jr  =  equivalent  damping  factor  for 
a  full  vibration  cycle. 


Figure  13.  Ice  force  plot  generated  by  numerical  inte¬ 
gration  starting  from  the  ice  crushing  strength  curve 
( test  40,  v  =  10  mmfs). 


second  mode  (60  Hz)  was  also  unstable,  indicating 
the  initiation  of  this  mode  at  high  ice  velocities.  In 
two  tests  the  2.6-Hz  mode  was  stable;  there  did  not 
appear  to  be  any  resonant  vibrations  in  this  fre¬ 
quency.  Vibrations  at  2.2-2. 8  Hz  are  due  to  res¬ 
onance  with  the  carriage.  These  comparisons  with 
theory  and  measurements  are  additional  evidence 
that  ice-induced  negative  damping  is  the  reason  for 
ice-induced  structural  vibrations  and  that  both  the 
ice  and  the  structure  affect  the  observed  ice  forces. 

BUCKLING  LOAD 

Tests  50  and  51  give  a  check  on  buckling  load 
theories.  Even  though  no  analytical  or  Finite  Ele¬ 
ment  Method  (FEM)  solutions  are  presented  for 
the  exact  geometry  of  a  circular  pile  in  an  infinite 
ice  field,  the  FEM  model  by  Sodhi20  is  close  enough 
to  the  problem  in  question.  For  compressive  stresses 
against  the  pile,  frictionless  boundary  conditions 
and  a  modulus  of  elasticity  of  3  GPa  yield  buckling 
loads  of  27  kN  or  greater.  The  measured  buckling 
loads  were  3.2  and  3.7  kN  with  1 02-  and  1 78-mm 
pile  diameters,  respectively. 

Ice  index  strength  tests  were  still  being  developed, 
and  therefore  the  measured  values  for  the  modulus 
of  elasticity  are  not  very  accurate.  Instead,  the 
modulus  of  elasticity  can  be  calculated  from  the 
measured  buckling  load.  The  curves  of  Sodhi  shows 
a  characteristic  length  of  about  350  mm,  which 
corresponds  to  a  modulus  of  elasticity  of  0.22  GPa. 

Buckling  developed  in  2-3  seconds.  Thus,  the 
loading  rate  is  high  enough  so  that  viscoelastic  yield¬ 
ing  can  not  be  the  sole  reason  for  the  difference. 

The  linear  theory  itself  must  be  at  least  part  of  the 
reason.  Kerr22  has  proved  (but  not  solved  the  real 
case)  that  the  bifurcation  point  in  the  equilibrium 
path  of  the  plate  on  the  clastic  foundation  under 
in-plane  loads  is  unstable.  This  means  that  imper¬ 
fection  sensitivity  can  significantly  reduce  buckling 
loads.  Meaningful  comparisons  between  the  meas¬ 
urements  and  the  theory  will  become  possible  after 
the  nonlinear  buckling  problem  is  solved. 

In  both  buckling  cases,  the  ice  failed  by  buckling 
only  when  the  velocities  were  below  5  mm/s;  above 
this  velocity  only  crushing  failure  occurred.  For  a 
given  ice  sheet  the  test  with  a  small  pile  did  not  pro¬ 
duce  any  buckling,  e.g.  test  49.  This  indicates  that 
buckling  and  crushing  loads  were  almost  identical 
at  the  low  velocities  in  tests  50  and  51.  A  slight 
change  in  conditions  was  enough  to  change  the 
mode  of  failure  from  buckling  to  crushing.  With 
increasing  velocity  the  effective  modulus  of  elas¬ 
ticity  of  ice  increased  while  viscoelastic  effects  were 


decreasing;  the  foundation  modulus  also  increased 
due  to  increasing  hydrodynamic  support  effects. 

These  measurements  and  the  analysis  of  the  re¬ 
sults  suggest  that  for  wide  full-scale  structures  with 
slowly  moving  ice,  buckling  can  occur  before  crush¬ 
ing.  This  could  be  important  for  design  ice  forces. 
However,  as  the  nonlinear  buckling  model  is  missing, 
it  is  not  possible  to  locate  the  boundary  between 
crushing  and  buckling. 

CONCLUSIONS 

Model  tests  conducted  at  the  CRREL  Ice  Engi¬ 
neering  Facility  during  December  1979  and  January 
1980  provided  valuable  data  for  understanding 
dynamic  ice-structure  interactions  during  crushing. 
The  test  arrangements  allow  scale-model  studies  for 
most  bottom-founded  offshore  structures.  In  this 
study,  only  constant-thickness  icc  was  used,  but 
modeling  is  possible  for  rafted  ice  and  pressure 
ridges  as  well. 

In  these  test  arrangements  the  weakest  link  is  the 
inadequate  stiffness  of  the  ice-pushing  system.  The 
dynamic  response  of  the  ice-pushing  system  still  in¬ 
terfered  at  low  velocities,  even  though  the  whole 
velocity  range  from  ductile  to  brittle  could  be  meas¬ 
ured. 

Experiences  with  urea-doped  ice  confirm  that 
sea-ice  grain  structures  and  crushing  failure  types 
are  similar.  The  crushing  strength  and  the  modulus 
of  elasticity  have  the  same  ratio.  Hence,  scale  model 
measurements  can  be  made  that  arc  accura’c  and 
simple  to  interpret. 

The  general  behavior  of  the  icc  force  and  pile 
vibrations  was  such  that  at  low  velocities,  a  sawtooth 
ice  force  and  pile  displacement  pattern  developed. 
With  increasing  velocity  the  lowest  natural  mode  of 
the  pile  was  excited;  it  dominated  the  response  until 
the  velocity  became  too  high  for  the  first  mode  to 
follow.  Then  cither  higher  frequency  ice  forces 
appeared  or  the  second  mode  of  the  structure  be¬ 
came  dominant.  The  shift  from  the  sawtooth  pattern 
to  the  natural  modes  of  the  structure  occurred  in 
each  test. 

Maximum  icc  forces  usually  occurred  at  low  ve¬ 
locities  in  the  range  of  loading  rale  from  ductile  to 
brittle.  The  maximum  and  average  ice  loice,  as  well 
as  the  standard  deviation,  tended  to  be  lower  with 
increasing  velocity,  although  the  random  peaks  of 
the  icc  force  at  high  velocities  may  equal  the  peaks 
at  low  velocities.  There  is  more  randomness  with 
freshwater  ice  than  with  urea-doped  ice.  I  or  design 
purposes  a  good  estimate  of  the  maximum  ice  torce 
is  the  average  plus  three  times  the  standard  deviation. 


These  tests  were  the  first  to  measure  the  effect 
of  the  ratio  of  pile  diameter  to  ice  thickness  on  the 
effective  ice  pressure  under  continuous  dynamic 
crushing.  The  results  indicate  that  a  similar  depend¬ 
ence  of  ice  force  on  aspect  ratio  under  static  condi¬ 
tions  holds  with  sufficient  accuracy.  It  appeared 
that  in  urea  ice  the  proportion  of  the  ice  layer  with 
a  horizontal  c-axis  more  strongly  affects  the  crush¬ 
ing  strength  than  the  bending  strength. 

The  frequency  of  the  observed  ice  force  is  directly 
proportional  to  the  ice  velocity  and  inversely  pro¬ 
portional  to  the  ice  force  fluctuation  range.  The 
natural  frequencies  of  the  pile  dominate  the  fre¬ 
quency  response  in  such  a  way  that  the  crushing 
frequency  stays  close  to  the  natural  frequency,  and 
instead  of  the  frequency  increasing  with  velocity,  it 
causes  the  ice  force  range  to  increase  with  velocity. 

The  test  results  support  the  theory  that  the  crush¬ 
ing  frequency  can  be  calculated  from  the  product  of 
the  structural  stiffness  at  the  ice  action  point  and 
the  ice  velocity  divided  by  the  observed  ice  force. 

The  dominating  feature  of  natural  frequencies  has 
to  be  taken  into  account  and  the  effective  structural 
stiffness  observed  if  applicable.  T  hus,  the  whole 
range  of  observed  frequencies  tan  be  described  by 
a  simple  equation. 

Theoretical  predictions  according  to  the  theory 
of  ice-induced,  self-excited  vibrations  yield  similar 
results.  Another  theory,  which  says  that  the  crush¬ 
ing  frequency  can  he  determined  by  the  property  of 
ice  to  fail  into  floes  of  a  certain  size,  did  not  agree 
with  the  test  data.  Peyton’s  suggestion  that  ice  has 
a  characteristic  failure  frequency  was  disproved. 

With  changing  velocity  the  crushing  frequency 
usually  changes;  it  remains  constant  only  at  the 
expense  of  the  ice  force  or  displacement  amplitudes 
and  is  the  result  ol  a  dominant  natural  mode  of  the 
structure. 

The  acceleration  data  reveal  that  during  crushing, 
high  frequency  components  were  present  and  gave 
the  highest  accelerations.  Altei  crushing,  the  struc¬ 
ture  settles  clown  to  first  or  second  mode  vibrations. 

While  the  ice  edge  sticks  to  the  pile  during  load 
build-up,  it  provides  additional  support,  which  in¬ 
creases  the  natural  liequencies.  Integrated  velocities 
verily  that  the  pile  deflection  >ate  is  so  high  that  the 
ductile  loading  rate  can  exist  dm  ing  the  whole  ve¬ 
locity  range,  despite  the  lad  that  the  ice  velocity 
alone  would  tesuli  in  brittle  loading  iates. 

Damping  measurements  indicated  that  ice  crush¬ 
ing,  friction  and  the  grinding  ol  ice  into  i  uhhle 
against  the  pile  cieated  considerable  positive  d. imp¬ 
ing.  mole  than  six  times  that  ol  the  pile  alone. 

However,  the  negative  damping  effects  ol  the  ice 
were  so  high  that  aperiodic  diver  gences  occurred 
at  the  beginning  ol  c  rushing.  As  this  was  only  dui  mg 
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a  part  of  the  vibration  cycle,  the  average  damping,  or 
the  equivalent  viscous  damping,  was  still  positive. 

The  origin  of  ice-induced  vibrations  was  verified. 
The  negative  damping  caused  dynamic  instability 
in  the  small,  but  as  the  overall  damping  was  positive, 
the  steady  state  vibrations  stabihzed  and  the  pile 
reached  its  limit  cycles.  Ice  force  histories  for  struc¬ 
tures  similar  to  those  measured  can  be  calculated 
using  the  ice-crushing  strength  vs  loading  rate  curves 
and  adjusting  for  the  dynamic  response.  This  gives  a 
basis  for  better  design  against  dynamic  ice  forces. 
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APPENDIX  A:  CORRECTED  ICE  FORCE  VS  TIME  FOR  TESTS  39  AND  44 
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